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We show that metamaterials can be used as a testing ground to investigate spontaneous symmetry breaking 
associated with non-Hermitian quantum systems. By exploring the interplay between near-field dipolar coupling 
and material absorption or gain, we demonstrate various spontaneous breaking processes of the PT -symmetry 
for a series of effective Hamiltonians associated to the scattering matrix. By tuning the coupling parameter, 
coherent perfect absorption, laser action and gain-induced complete reflection (n reflector) by using an ultra- 
thin metamaterial can be obtained. Moreover, an ideal PT -symmetry can be established effectively in a passive 
system by using metamaterials. 

PACS numbers: 78.67.Pt,42.25.Bs,78.20.-e 


Metamaterials, periodic arrangements of artificial subwave¬ 
length atoms, have paved the way for realizing different novel 
concepts and applications in optics, such as negative re¬ 
fraction, sub-diffraction imaging and invisibility cloaking|l- 
[5]. Instead of using the macroscopic material parameters 
for description, many of the optical properties, including 


metamaterial-based chirality 
duced transparency (EIT)l 


and electromagnetic ally in- 
, can be understood by study¬ 


ing the coupling between different magnetic and electric artifi¬ 
cial atoms in the so called hybridization picture Bisr 
Similar to optical lattices with controllable coupling | 
the design flexibility in the coupling between artificial atoms 
opens up a new perspective to realize quantum simulators 
to achieve effects which are difficult to obtain in quantum 
systems fl7l fljll . Here, we theoretically investigate the pos¬ 
sibility in constructing a metamaterial system with an effec¬ 
tive PT -symmetric Hamiltonian which can still show real 
eigenenergies. One of the most interesting properties of a 
PT -symmetric Hamilton is the existence of a phase tran¬ 
sition arising from spontaneous symmetry breaking by tun¬ 
ing an external parameter, defined as the coalescence of two 
eigenvalues and their corresponding eigenvectors |(20i]. Un¬ 
like the previously considered symmetry breaking effects of 
metamaterials §21], such a spontaneous symmetry breaking is 
implicit in nature. It means that the phase transition occurs 
without structurally destroying a specific symmetry (e.g. mir¬ 
ror, rotational symmetry). The external parameter being tuned 
still respect the PT -symmetry. 


During the past few years, PT -symmetric Hamiltoni¬ 
ans have in fact been extensively investigated within the 
framework of optics. The most direct analog is repre¬ 
sented by a pair of optical waveguides with a balanced gain 
and loss profile 1191 122143 ill . Spontaneous PT -symmetry 
breaking has also been theoretically and experimentally 
investigated0 I22U3H1 . Intriguing phenomena, such as 
asymmetric transmission !!^!. El , power oscillation [El, uni¬ 
directional invisibility ll 1 9L 2911 . coherent perfect absorption 
and laser[30] have been predicted. Moreover, extending from 
an ideal PT -symmetric system with balanced gain and loss, 
non-ideal PT -phase transition of a passive system has also 


been proposed and experimentally demonstrated E 0. 
Unconventional loss-induced optical transparency is demon¬ 
strated around the phase transition point[25]. Actually, PT - 
symmetry has not been discussed within the context of meta¬ 
materials. By employing metamaterials, the larger flexibility 
in designing the coupling between different artificial atoms 
gives us additional possibilities in realizing non-trivial PT - 
symmetry. Here, we show that by mapping a metamaterial 
consisting of a bright electric and a dark magnetic atom to 
an effective Hamiltonian, we can establish an ideal sponta¬ 
neous symmetry breaking transiting from little absorption to 
complete absorption even in a passive system. It relies on the 
inverse relationship that scattering is much larger than loss for 
a bright electric atom while the reverse is true for a dark atom. 
Such a system has been previously considered to achieve EIT. 
However, we are interested in the phase transition to obtain 
complete absorption[32l El across a threshold in the cou¬ 
pling parameter, in contrast to the high transmission for EIT 
emerging from the strong coupling regime |0|- Moreover, by 
incorporating a small material gain into the dark atom of the 
ultrathin metamaterial, we can obtain laser action lE -Ebll and 
gain-assisted complete reflection by considering the sponta¬ 
neous symmetry breaking effect of a class of Hamiltonians 
associated with the scattering matrix. 

The discussion of PT -symmetric Hamiltonians originates 
from Bender and colleague’s work, showing that a Hamilto¬ 
nian H respecting a parity-time symmetry, i.e. [PT, H] = 0, 
can have completely real eigenenergies [s2Q]. To find the ana¬ 
log of PT -symmetry in a metamaterial system, the simplest 
way is to look for a pair of coupled artificial atoms, since 
their response can be understood within a hybridization pic¬ 
ture in quantum mechanics as a two state problem, and the 
PT -symmetry in two-state quantum system is already well 
studied 123 ■ Here, we consider a metamaterial constituting a 
planar array of atoms. Each unit cell consists of an electric 
atom coupled to a magnetic atom through the near fields. For 
simplicity, we assume normal incidence and there is no polar¬ 
ization conversion for the incoming plane waves. The optical 
response of such a metamaterial can be expressed using the 
scattering matrix S as 
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where a± and b± are the E-field amplitudes of the incoming 
and outgoing waves in the positive and negative directions. 
The same scattering process has a microscopic representation 
using the atomic responses and their associated far-field ra¬ 
diations. By further assuming each atom is dominated by 
only one dipolar resonating mode, the atomic responses can 
be summarized as 


(Ho - col) 


P 

im 


^y\(a + +aJ) \ 

i yfyl( a + ~ a -) / 


(2) 


while the outgoing fields and the atomic responses are related 
to each other by 


(H 2 - ml) 


P 

im 
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In the above equations, y,- = y s + y ,oss (i = 1,2) represents 
the total resonating linewidth for each atom where y s , a pos¬ 
itive number, represents the part of scattering and y ,oss is the 
part due to absorption loss when it is positive. We also nor¬ 
malize the electric dipole p and the magnetic dipole iii by the 
oscillator strengths as in Ref.7, and assume the two resonant 
frequencies at the same frequency coi res — too- In fact, Eq. (4) 
defines a series of Hamiltonian H n with meaning in the fol¬ 
lowing. By taking an appropriate arithmetic average of Eq.(2) 
and Eq.(3), one can prove that det(H n - ml) = 0 yields an 
eigenvalue (n - 2 )/n of matrix S. We can also obtain 


det(H 2 - ml) _ det(Hi -ml) 
det(Ho - ml) det(Ho - ml) 


(5) 


For example, det(Ho - ml) = 0 leads to the pole of S, which 
corresponds to laser actionlflol l34l - [36ll . det(H 2 - ml) = 0 
implies a zero of S to completely absorb the incoming EM 
waves, namely coherent perfect absorption (CPA)B0J32ll331. 
For a metamaterial with mirror symmetry in the propagation 
direction, the two eigenvectors of S are (a 4 .,a ) = (1,1) with 
eigenvalue t+r and (1, -1) with eigenvalue t—r by defining t/r 
as the complex transmission/reflection coefficient for electric 
field with respect to the mirror plane. If we further confine our 
discussion to the coupling between a bright electric and a dark 
magnetic atom, we have the additional requirement t = 1 + r. 
Then, det(Hj - ml) = 0 implies a complete reflection r — — 1 
and t = 0. In general, a particular n/eigenvalue of S denotes 
a particular optical phenomenon and it corresponds to solv¬ 
ing the eigenvalue problem of H n . If the eigenvalue of H n 
is real, the phenomenon can occur at a real frequency with a 
pole/zero on the corresponding optical spectrum. Therefore, 


by considering the spontaneous symmetry breaking process 
of a PT -symmetric H n , one can observe the phase transi¬ 
tion when the real eigenfrequencies are turned into complex 
conjugates or vice versa. For our 2x2 Hamiltonian, the P- 
operation switches rows and columns of the matrix while the 
T-operation takes complex conjugate on the elements. A PT- 
symmetric H n is defined by [PT, H n ] = 0. We have to empha¬ 
size that PT -symmetry being discussed here is an effective 
one through the establishment of an effective PT -symmetric 
Hamiltonian, which is different to the approach employing 
complex potentials |20h . 

As an example, a PT -symmetric H 2 corresponds to the 
condition y 2 - 2y s 2 = 2y\ - y\, which can be satisfied by a 
representative hybridized metamaterial with a bright electric 
atom (scatters most) and a dark magnetic atom(absorbs most). 
For a fixed set of parameters, e.g. y\ = 15 THz, y l ° ss = 1.75 
THz (yj = 13.25 THz), y 2 ~ y'° ss , this condition is repre¬ 
sented by the vertical red line in Fig. 1 at y 2 = 11.5 THz 
to indicate a PT -symmetric H 2 . Then, the spontaneous sym¬ 
metry breaking process can be observed by changing the cou¬ 
pling coefficient k (still respecting the PT -symmetry) across a 
threshold value. The eigenfrequencies of PT -symmetric H 2 
are reaKsolid red line) when k > y 2 , and becomes a com¬ 
plex conjugate pair(dotted red line) when k < y 2 . The split¬ 
ting of the real/imaginary part of the eigenfrequencies with 
varying k is showed in the inset. This behavior reveals itself 
through the absorption coefficient when the metamaterial is 
excited by the incident waves of (a+, a_) = (1,1), giving CPA 
for k > y 2 at the two real frequencies. We note that such 
a configuration of coupling between a bright and a dark atom 
(without the PT -symmetric condition satisfied in general) has 
been used to give EIT in the strong coupling regime 01. We 
endow the same system to visualize effective PT -symmetry 
breaking and focus on its optical properties in the regime near 
the critical point. In the following, we demonstrate the inter¬ 
esting optical properties associated with PT -symmetry break¬ 
ing with metamaterial structures. In particular, undergoing the 
phase transition, the metamaterial becomes perfectly absorp¬ 
tive with real frequencies. The bright atom (inset of Fig. 2(a)) 
is a silver strip of length LI = 150 nm, width W 1 = 50 nm 
and thickness t — 20 nm. Fig. 2 (a) shows the electric field 
near one end of the strip with numerical simulation (symbols) 
and theoretical model (lines). The numerical simulation is car¬ 
ried out using finite element analysis (Comsol MultiPhysics) 
with Drude model of silver being e - - or p (u) 2 + iy p cS)~ x , 

where io p = 2.1961 x 10 15 Hz, y p = 4.3439 x 10 12 Hz, and 
too = 3.7. The theoretical result is fitted using a Lorentz model 
E s oc (coo -co- fyi) 1 (yi - y l ° ss )Eo, where a>o = 402.5 THz, 
yi = 15 THz, and y'° ss = 1.75 THz, E m is the scatter- 
ing(incident) E field. y‘° ss can be estimated by the probing the 
response of the bright atom with and without metal loss. The 
dark atom is a magnetic resonator tuned to the same resonat¬ 
ing frequency with two parallel silver strips and a dielectric 
gap material (e = 2.25 + ie" with e" introduced to represent a 
slight absorption or gain) in between. The geometric param¬ 
eters are L2 = 100 nm, W2 = 30 nm, f = 20 nm, as showed 
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FIG. 1: (Color online) Parametric evolution of the metamaterial for 
different H„(n = 0,1,2); Inset: Splitting of eigenvalues for H 2 with 
varying k. Vertical solid lines indicate the real eigenvalues case for 
ideal PT -symmetric H„ above the phase transition threshold, short 
dotted lines illustrate the complex conjugate pair eigenvalues case for 
ideal PT -symmetric H n below the phase transition threshold. Inset 
gives the eigenvalues with varying k for PT -symmetric H 2 . Curved 
solid lines show the real eigenvalue condition indicating the optical 
properties for non-ideal PT -symmetric H„, short dashed lines show 
branching condition for non-ideal PT -symmetric H„. 


in the inset of Fig. 2(b). Due to the large quality factor of 
the dark mode resonance, a small e" of the gap material can 
effectively alter the half-linewidth y 2 . As shown in Fig. 2 (b), 
numerically calculated field (symbols) near one of the metal 
strip ends is again fitted well with theoretical Lorentz model 
(solid line) for e" = 0,-0.1 (hollow, solid symbols). Com¬ 
paring between numerical and analytical results with different 
e ", a simple approximated formula between yi and e" can be 
established as 72 = 1-75 + 22.5e". Now, we hybridize these 
two atoms within a metamaterial with a square lattice of pe¬ 
riodicity d = 400 nm and a distance h in separating the two 
atoms, as shown in Fig. 2(c). The change of distance h be¬ 
tween the two meta-atoms can tune the coupling k between 
them, but there exists no simple relationship between them, as 
shown in Fig. 2(d). The k can be calculated from the bright 
meta-atom probe response for e" - 0 by changing the dis¬ 
tance h between them. The fitted formula can represent main 
variation character of k, but exhibits slight deviation in large k 
region. 

Now, we illustrate the spontaneous symmetry breaking with 
fixed geometry of the meta-atoms with constant a> 0 , y 1 , and 
yioss. | | le variable parameters are the coupling coefficient k 
(or equivalently h) between the two atoms and 72 (through 
the tuning of e"). We can then choose a particular 72 to sat¬ 
isfy PT -symmetry requirement of H n . For a PT -symmetric 
H 2 , we need 72 = 2y\ - 71 = 11.5 THz, i.e. e" = 0.433, 
which is a passive configuration to display the effective PT 
symmetry. In the light of the above theoretical discussion, the 



FIG. 2: (Color online) Analytical (solid line) and numerical (scatter 
line) spectral response of the E-field probe for bright(a) and dark(b) 
meta-atom to establish the atomic Lorentz model with fitted param¬ 
eters (oj 0 THz, 7 ] THz, 7 |""', and y 2 ). (c) and (d): The hybridized 
metamaterial under consideration and the relationship between the 
coupling coefficient k and the distance li between the two meta¬ 
atoms. 



FIG. 3: (Color online) Observation of effective spontaneous PT - 
symmetry breaking for H 2 and Hi. Analytical model (with fitted 
parameters from Fig. 2) (a) and numerical(b) outgoing |t + r| 2 spectra 
versus frequency and coupling distance. Analytical(c) and numeri- 
cal(d) transmission spectra versus frequency and coupling distance. 

eigenvalue problem for a PT -symmetric H 2 is equivalent to 
the zeros of the S, which can be easily observed in the co¬ 
herent outgoing |r + t\ 2 spectra when we change the param¬ 
eter h. The PT -phase transition point in this case occurs at 
Kpj = 11.5 THz (h = 46.0 nm). Below the critical point, 
Kpj < 11.5 THz (h > 46.0 nm), the eigenvalues of H 2 is 
a complex conjugate pair, there are no crossing in real fre¬ 
quency axis, so we cannot observe CPA in this region, shown 
in Fig. 3(a). As we increase (decrease) /r(/z), at Kpp = 11.5 
THz(/z = 46.0 nm), the two eigenvalues of H 2 collapse to 
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one point in the real frequency axis, i.e. CPA(|r + t\ = 0) 
emerges in the spectra. When k passes over the Kpj, there 
are two real frequencies satisfies the eigenvalue solution of 
H 2 , CPA branching phenomenon can be clearly shown in Fig. 
3(a) from the analytical results obtained from Eq. (2) and (3) 
with the atomic model established from the previous simula¬ 
tions shown in Fig. 2. The same phenomenon can be ob¬ 
served from the numerical simulations (with microstructures 
included) shown in Fig. 3(b). Real parts (white dashed line in 
Fig. 3(a),(b)) of the eigenvalues of H 2 agree well with the CPA 
positions in the spectra. In both Fig. 3(a) and Fig. 3(b), the 
absorption loss (1 - \t + r\ 2 for symmetric incoming waves) of 
the hybridized system increases dramatically with increasing 
coupling (decreasing separation) below Kpj. At Kpj, the pres¬ 
ence of m interferes with p to the right size and gives r — — 4 
and t — \. This is the onset of CPA phase. Going beyond 
the transition point (larger coupling), the CPA peak splits into 
two. On the other hand, this splitting makes both atoms can¬ 
not be efficiently excited at u>o and loss decreases at that fre¬ 
quency. In fact, the non-monotonic function of loss against 
coupling coefficient is a direct manifestation of spontaneous 
symmetry breaking, which is similar to the loss-induced opti¬ 
cal transparency in a passive VT -symmetric system deviated 
from the ideal VT -symmetry [25]. Here, using metamaterials, 
we can actually establish an ideal VT -symmetry in a com¬ 
pletely passive system. CPA becomes possible besides the 
counterintuitive trend of loss function. 

Besides H2 in revealing VT -symmetry breaking process. 
Hi is also a potential candidate to display the same process. 
As predicted in Eq. (5) and Fig. 1, by taking into account 
the eigenvalue problem of Hi, we need to fulfill the case 
72 = -y,"” = -1.75 THz (the black vertical line in Fig. 1), 
i.e. e" = -0.156, where potential amplifying media can be 
available in nature The emergence of effective VT - 

symmetry breaking process can be observed though transmis¬ 
sion spectra, and t = 0 (r = -l)in spectra implies real fre¬ 
quency solution for eigenvalue of Hj. It indicates the system 
without absorption loss, A — 1 - \t + r\ 2 = 0, and transmis¬ 
sion, t — 0 , with perfect reflection with and a /r-phase shift, 
r = — 1. The phase transition point in this configuration is at 
Kpj- = 1.75 THz (h = 98.4 nm). Illustrated in Fig. 3(c), and 
numerically verified in Fig. 3(d), the phase transition point 
is at h = 97.5nm, t = 0 (r = - l)in the transmission spec¬ 
tra agrees well with the real parts of the eigenvalues of Hi. 
Below the phase transition point (k < Kpj), the eigenvalue of 
H\ implies no cross point in the real frequency axis, a com¬ 
plex conjugate pair, so transmission is dominated by the bright 
atom dipole oscillation, although t as 0, the eigenvalue of S is 
not -1 (Absorption loss is not zero). Above the phase transi¬ 
tion point Kpj, the eigenvalue of H\ implies two real eigen- 
frequencies, sof = 0(r = -l) can be observed at two differ¬ 
ent frequencies in the transmission spectra. Real parts (white 
dashed line in Fig. 3(c),(d)) of the eigenvalues of Hi match 
the t = 0 positions well in spectra. The effective spontaneous 
VT -symmetry breaking for Hi can thus be observed in the 
spectra. It is worth to note that the transmission exhibits maxi¬ 


mum value at the position k s = y -7271 = 5.12 THz (h = 73.8 
nm for analytical prediction and h — 72.5 nm for numerical 
results) at a>o, because gain behavior of the system appears, 
when pass over the phase transition point Kpj-- As we in¬ 
crease k , we can encounter a real eigenvalue of Ho , indicating 
pole of S, which implies the_system shows laser action at this 
parameter position DOl 134143611 . as indicted in the cross point 
between the black vertical line indicating VT -symmetric Hi 
and the blue curve displaying real eigenvalue of Hq in Fig. 1. 
The laser action results from non-ideal VT -symmetry of Ho 
whose origin will be discussed in the following. 

For Ho, ideal VT -symmetry condition 72 = -71 implies 
that huge gain is required to compensate metal losses. This 
is difficult to be satisfied using natural materials. However, 
in a practical (non-ideal) system, even though ideal condi¬ 
tion is not satisfied, the eigenvalues still split in the same way 
as an ideal VT -symmetric Hamiltonian. Only an additional 
constant bias in the imaginary part of the eigenvalues occurs. 
Similar strategies have been employed in studying non-indeal 
VT -symmetry in passive systems H25L 27]. Then, phase tran¬ 
sition (branching) point is not on the real frequency axis, and 
is shifted into the lower-half of the complex w-plane. One of 
the eigenvalues (the one with increasing imaginary part) will 
cross the real frequency axis, corresponding to the original 
optical phenomenon, e.g. laser action for the case of Ho as a 
pole of S at k s = y— y 1 y 2 130ll34l - i36[| . It is the solid blue curve 
in Fig. 1 and is lying below branching condition (branching 
occurs when k > Kpj) where Kpj = (71 - 72 )/ 2 (dashed green 
line). In contrary, ideal VT -symmetric Ho has the poles of S 
on the real frequency axis beyond branching condition (ver¬ 
tical solid blue line) at two different frequencies. As an ex¬ 
ample, 72 = -0.5 THz (e" = -0.1), Fig. 4 gives the outgo¬ 
ing intensity of the system with symmetric incoming light as 
a function of k(/i). Laser action collapses to one position at 
a> 0 , h = 89.4 nm for theoretical prediction and h = 90.0 nm 
for numerical results. The branching phenomenon still occurs 
(without real poles of S) and agrees well with the real parts 
of the eigenvalues of Ho (black dashed lines), as displayed 
in Fig. 4. Laser action and branch phenomenon is the direct 
demonstration of non-ideal 'P7~-symmetry for Ho. Besides 
Ho, non-ideal VT -symmetry can also occur for Hi and H 2 . 
In Fig. 1, the real eigenvalue of H 2 (zero of S), i.e. CPA, 

is at k 0 = ^ 72(71 - 2y ! ° ss ) (solid red curved line) while the 
branching condition is at Kpj = (71 + 72 - y l ° ss )/2 (dashed 
purple line). For Hi, perfect reflection, r = —1, t = 0 (one S 

eigenvalue is - 1 ) is at k^\ = ^j-j 2 'y , ° ss (solid black curved 
line) while the branching condition is at Kpj = ( y l ° ss - 72)/2 
(dashed black line). So, for a non-ideal VT -symmetric sys¬ 
tem, the optical phenomenon (e.g. laser, CPA or ^--reflector) 
associated to H„ collapses to a single point (with a k a hit less 
than the one triggers branching) instead of the two branches. 

In the above discussion, k is the external parameter in 
demonstrating VT -symmetry. If we change 72 with a fixed 
k instead (travelling in a horizontal direction in Fig.l), the 
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h (rim) h (nm) 

FIG. 4: (Color online) Output intensity (color map) of the metama¬ 
terial log w \t + r\ 2 revealing spontaneous symmetry breaking of non¬ 
ideal VT -symmetry for Ho and the real part of the eigenfrequencies 
for Ho(dashed lines), (a) Analytical model; (b)Numerical simula¬ 
tions with metamaterial microstructures. 



FIG. 5: (Color online) Analytical (solid line) and numerical (dashed 
line) outgoing spectra, as a function of y 2 for the metamaterial at oj 0 . 


metamaterial can function as a laser, the cross point for Ho, 
or a CPA, the cross point for H 2 at ojq. As an example for 
demonstration, h = 75 nm, besides appearance of laser ac¬ 
tion and CPA in spectra, schematically shown in Fig. 5, above 
CPA position, an increase of loss (increasing 72 ) leads to more 
output of the system, which is similar to loss induced optical 
transparency in a noil-ideal passive TT system ll25tl . Similar 
counter-intuitive trend also occurs beyond the laser position, 
an increase in gain (decreasing negative yi) causes less output 
observed in the spectra. 

In conclusion, we investigated the effective VT -symmetric 
Hamiltonians using a metamaterial system to reveal sponta¬ 
neous symmetry breaking processes. Theoretical and numer¬ 
ical results clearly confirm the analog of spontaneous break¬ 
ing of VT -symmetry using metamaterials. Our results not 
only open up a new perspective towards testing the intrigu¬ 
ing properties of spontaneous VT -symmetry breaking in non- 
Hermitian quantum systems, but also provide feasibility for 
future applications in plasmonic and metamaterial systems, 
such as laser action i30ll34l43611 . and CPA with sub-wavelength 
thickness [33- Moreover, by using metamaterials, it becomes 
possible to establish an ideal VT -symmetry by using a passive 
system. 
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